非齐次动力学方程的一种精细积分单步方法

针对非齐次动力学方程■,结合精细积分法和微分求积法,利用同阶的显式龙格-库塔法对计算过程中待求的v_(k+i/s)(i=1,2,…,s)进行预估,提出了一种避免状态矩阵求逆的高效精细积分单步方法。该方法采用精细积分法计算e~(Ht),而Duhamel积分项采用s级s阶的时域微分求积法,计算格式统一且易于编程,可灵活实现变阶变步长。仿真结果表明,与其他单步法及预估校正-辛时间子域法进行数值比较,该方法具有高精度、高效率及良好的稳定性,在求解大规模动力系统时间响应问题中具有较大的优势。     

形状记忆合金椭圆孔口问题的有限单元法

基于Lagoudas形状记忆合金(SMA)三维本构模型,假设材料为各向同性,推导了SMA平面应力状态的增量型本构方程,继而编写了ABAQUS用户自定义材料(UMAT)子程序,研究了在双向拉伸情况下,外载荷、温度、椭圆孔口长短轴之比对超弹性SMA椭圆孔口板中应力诱发马氏体相变区的影响。数值结果表明:应力诱发马氏体相变首先发生在椭圆孔口长轴端点部位,在外加载荷作用下逐渐扩展到板内,并由内向外形成马氏体相区、相变混合区和奥氏体相区;SMA板内应力诱发马氏体完全相变区面积与施加外载荷成正相关,与温度成负相关;随着椭圆孔口长短轴之比增大,SMA板内应力诱发马氏体完全相变区面积呈现出先减小后增大的趋势;拉应力差值相同时,相较于拉应力沿椭圆孔口长轴方向较大的情况,当拉应力沿椭圆孔口短轴方向较大时,SMA板内完全相变区面积较大,椭圆孔口周边应力集中现象更明显。     

基于图像的二维剪纸自动生成方法

中国剪纸的设计极具挑战性, 要求画面简洁、直观, 还需要表达特定的文化内涵, 且整张剪纸须整体连通。提出了一种基于图像的二维剪纸自动生成方法, 能够将任意数码照片自动转化为剪纸图形。首先，利用图像分割方法建立区域连接图; 接着, 基于该连接图对颜色、边界对比度和区域连通性进行数学建模, 并获得优化目标函数; 最后, 通过模拟退火算法求解目标方程, 自动生成保持图像内容的剪纸图形。还开发了连通性后处理和区域指定等用户交互工具, 允许用户在自动生成的剪纸图形中方便地加入个人设计。实验表明, 所生成的剪纸图形画面简洁、整体连通。 本方法在降低剪纸设计难度的同时还可满足个性化的设计需求, 有助于传播和传承我国的民间剪纸艺术。     

介质材料二次电子发射特性对微波击穿的影响

以介质填充的平行板放电结构为例,本文主要研究了介质填充后微波低气压放电和微放电的物理过程.为了探究介质材料特性对微波低气压放电和微放电阈值的影响,本文采用自主研发的二次电子发射特性测量装置,测量了7种常见介质材料的二次电子发射系数和二次电子能谱.依据二次电子发射过程中介质表面正带电的稳定条件,计算了介质材料稳态表面电位与二次电子发射系数以及能谱参数的关系.在放电结构中引入与表面电位相应的等效直流电场后,依据电子扩散模型和微放电中电子谐振条件,分别探讨了介质表面稳态表面电位的大小对微波低气压放电和微放电阈值的影响.结果表明,介质材料的二次电子发射系数以及能谱参数越大,介质材料的稳态表面电位也越大,对应的微波低气压放电和微放电阈值也越大.所得结论对于填充介质的选择有一定的理论指导价值.     

智能蜂窝状有序多孔薄膜：体系构建、响应性能及应用探索

蜂窝状有序多孔膜是一种典型的具有微、纳米尺寸的图案化高分子薄膜材料,在超疏水表面、光电材料、组织工程、生物医学、传感器等领域都具有良好的应用前景。表面形貌或性质具有可逆刺激响应特性的智能型多孔膜不仅拥有常规多孔膜的性能,而且还拥有独特的"开-关"功能,因此受到广泛关注。本文回顾了智能有序多孔膜近年的研究进展,重点对其体系构建、响应方式、机理、性能和相关应用进行了归纳与分析,评述了不同刺激响应多孔膜的特点及优劣势,探讨了智能多孔膜的优化设计及提高响应性能等问题,指出了目前面临的挑战,并展望了未来的发展方向。     

浸没边界-简化热格子Boltzmann方法研究及其应用

针对流固耦合传热问题,本文提出了一种基于浸没边界-简化热格子玻尔兹曼方法(immersed boundary method-simplified thermal lattice Boltzmann method,IB-STLBM)的耦合模型.不同于传统的格子玻尔兹曼方法使用分布函数演化流场和温度场,简化热格子玻尔兹曼方法(simplified thermal lattice Boltzmann method,STLBM)的演化过程不需要依赖分布函数,只涉及平衡态分布函数和非平衡态分布函数,能够直接演化宏观量,极大减小了计算过程中所占用的虚拟内存,简化了边界条件的实现方式,同时具有较高的稳定性.传统的浸没边界法对流场的计算采用欧拉网格,对固体边界采用拉格朗日网格,认为固体边界是对流场产生某种体积力.在应用浸没边界法时,汲取介观的思想,把固体的介入看作是对流场的干扰,打破了固体附近流体介观微团颗粒原始的平衡状态,这种干扰可以看作是在耦合边界上产生的一个非平衡项,可用非平衡态分布函数来表示.基于此,在模型中浸没边界法与简化热格子玻尔兹曼方法更紧密联系在一起,更大程度发挥二者的优点,整个计算过程更加简单直观,符合物理特性.通过对热圆柱绕流和内含热颗粒的封闭方腔自然对流问题的模拟以及对其结果的分析,验证了该算法在求解流固耦合传热问题的有效性和可行性.      

砂土类材料离散元细观参数的自动识别

准确的细观参数是采用离散元进行合理模拟的前提,基于最优化理论,本文提出了一种砂土类材料离散元细观参数的自动标定方法:采用PFC离散元软件模拟砂土的双轴压缩试验,以使PFC数值模拟的应力–应变曲线与某标准砂土室内双轴压缩实验结果的误差最小为优化目标,采用模式搜索法实现砂土类材料离散元细观参数的自动标定。本文提出的离散元细观参数自动标定方法具有识别效率快,精度高的优点。     

基于DIC的爆炸加载下脆性材料裂纹扩展规律的试验研究

通过数字图像相关法（DIC），应用PMMA对爆炸加载条件下脆性材料的裂纹扩展规律进行了试验研究。基于对称性试验模型，实现了裂纹尖端位置和应变场信息的同步记录。以此为基础，通过对比分析获知，主应变场应变值最大点不能作为裂纹尖端的判断依据。并以动态裂纹扩展速度为参量，应用断裂动力学和最小二乘牛顿迭代法，计算出了考虑惯性效应的Ⅰ-Ⅱ混合型裂纹的应力强度因子:K_Ⅰ和K_Ⅱ值会随着裂纹扩展方向改变而发生突变；K_Ⅰ最大值为2.63 MPa·m1/2，最小值为0.89 MPa·m1/2；其整体变化趋势表明，爆炸加载条件下脆性材料裂纹扩展随能量积聚和释放呈循环阶梯式递减发展。     

考虑空化效应的螺旋槽机械密封液膜动力学特性研究

考虑液膜空化效应的影响，研究螺旋槽液体润滑机械密封的动力学特性. 基于液体润滑理论和小扰动法，建立了考虑液膜空化的密封微扰膜压控制方程，采用有限单元法对端面液膜三自由度微扰下的液膜刚度和阻尼系数进行了数值求解，分析了不同参数对液膜密封动力系数的影响. 螺旋槽深度在10 μm左右、槽坝比在0.75左右、槽宽比在0.4左右，螺旋角在9°左右时液膜具有最大的轴向和角向刚度系数. 螺旋槽深度在5 μm左右、槽宽比在0.6左右、螺旋角在20°左右时，两角向交叉阻尼绝对值最大. 初始偏角的存在使密封压力呈现非对称性，从而使两角向动力系数绝对值不再相等. 液膜轴向刚度k_zz在数量级上远大于其余液膜刚度值，液膜轴向阻尼d_zz、角向阻尼d_αα和d_ββ远大于液膜其余阻尼值且随着转速和间隙的增大而减小.